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A different approach to the Fresnel region field approximation is introduced. Instead 
of using the conventional truncated power series expansion approximation, the Newton's 
iteration formula for square root is used. By using such an approximation for circular 
apertures with tapered illumination of the form (1 — £* 2 ) n , the Fresnel region fields can finally 
be expressed in terms of a new class of functions T^S(y> u )- The function Wq is shown to 
be related to the Lommel's functions of two variables, and the function Wl is then obtained 
from the function Wq~ 1 bv a simple recurrence relationship. Field distributions for 
n = 0, 1, 2, 3, 4 and at distances }{ D 2 /\, }i D 2 ;\, y 2 D 2 /\, D 2 /\, 2D 2 /\, oo have been computed 
and presented as sets of curves. General and quantitative properties of the fields are 
clearly demonstrated by these curves. It is also shown that the field of any other nonuniform 
illumination with circular symmetry can be expressed in terms of the fields of the basic 
illumination of the form (1 — f 2 ) n . 

1. Introduction 

It is well known that there is a close similarity between diffraction problems in optics and 
aperture field problems in antenna theory. In optics, two diffraction problems — the diffract ion 
from a single edge or a slit and the diffraction from a circular aperture — have been well treated. 
These two problems were first solved by Fresnel and Lommel, respectively, in terms of the 
now so-called Fresnel integrals and Lommel's functions of two variables. The corresponding 
problems in antenna theory are the Fresnel region fields produced by rectangular and circular 
apertures. However, in optics, the primary illumination is assumed to be uniform over the 
aperture, but in antenna theory, the case of nonuniform illumination is of great interest. 

For the case of uniform illumination, the analytical methods developed in optics can also 
be applied to antenna problems. For Fresnel region fields of rectangular apertures, generali- 
zations to include nonuniform illumination through expansions in terms of different basic 
forms are not too difficult, polynomial and trigonometric sum types of illumination were 
treated by the author [l]. 2 On the other hand, the analysis of Fresnel region fields of circular 
aperture antennas with nonuniform illumination is considerably more difficult. However, 
with the use of illumination of the form (1 — f 2 ) n , the Fresnel region field was solved by the 
present author [1]. A short communication about this result has been published [2]. The 
complete analysis and numerical results are presented in this paper. A refined but more 
complicated result for the case of uniform illumination was given by Hansen and Bailin [3]. 
It may be noted that the corresponding far field problems have been treated by many authors 
with the use of different basic forms [4, 5, 6, 7], including the form (1 — $ 2 ) n [4, 5] used in this 
paper. 

2. Fresnel Region Approximation 

In the study of field distributions of linearly polarized uniform phase plane aperture 
antennas [8] the following formula [1] or a slightly different form given by Silver [8] may be 
used as a starting point. 

E( w )^J A FfM,^(l+^dA, (l) 



1 Contribut Ion from Department of Electrical Engineering, Syracuse University, Syracuse 10, N.Y. This work is part of a project sponsored 
by Rome Air Development Center under Contract No. AF 30(f>02)-928 with Syracuse Univ. 
-' Figures In brackets Indicate I he literature references at the end of this paper. 
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where 



E(x,y,z) = the complex magnitude (amplitude and phase) of the electric field intensity at the 
field point (x,y,z), 
F(%,ri)=ihe illumination at the point (£,i?,0) in the apperture A. For uniform phase illumi- 
nation, F(^j7i) can be taken as a real function of £ and 77, 
r =^'(x— £) 2J r(y— 77) 2 +;2 2 =the distance between the field point (x,y,z) and the source 
point (£,17,0), and 






with X=the wavelength. 



The coordinate system used for writing formula (1) is given in figure 1. Formula (1) holds 
for any illumination function F(^,rf) and for apertures of any shape. In the present paper, it 
is assumed that the aperture is circular and the illumination is of circular symmetry. 

For a circular aperture of radius a, as given in figure 2, formula (1) can be rewritten as 
follows: 



where 



E(r,0, 0)=g fi'^Fip,*') ^p (l+£ cos ») pd P d<t>' 



r =-\J(r sin 6 cos <£— p cos </>') 2 + (? sin 6 sin <j>— p sin <t>') 2 + (r cos 0) 2 . 



(2) 



The evaluation of the above integral is a very difficult one without further approximation. In 
order to facilitate its evaluation, different approximations are generally used. Depending upon 
the approximation used, the result may be classified as the far field approximation, the Fresnel 
region approximation, etc. The most widely used approach to such approximations is to expand 
the distance r , appearing in the exponential, into a power series involving r, 6, and <j>. The 
result obtained, by neglecting all terms of order higher than one or of order higher than two, is 
generally considered as the far field approximation or the Fresnel region approximation. As 
far as the remaining r's under the integral sign are concerned, much more crude approximations 
are satisfactory. If the field point in consideration is not too close to the aperture and 6 is not 
too large, then the factor l/r (l + (r/r ) cos 0) can be approximated to good accuracy by using 
the relation, 

2 



-(l+-cos0W 



(3) 



U,y,z) 




Figure 1. Coordinate system for A plane aperture. 



(r, 0,</> ) 




Figure 2. Coordinate system for A circular aperture. 
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A rather different approach to the Fresnel region approximation is used in the present 
paper. Instead of using the truncated power series expansion, the well-known Newton's 
iteration formula for finding square root of a given number is employed. This formula has 
the following special property. If an approximation B lf of -yJA is known to be correct to n 
significant figures, then the approximation B 2 of -\JA obtained by using the iteration formula, 



a 



-K b ' + t) <« 



will be correct to 2n significant figures. If this formula is applied to the approximation of ro 
in Fresnel region, we have 

If (r sin 6 cos </>— p cos </>') 2 +( rsm # sm — P sin<//) 2 +(> cos 0) 2 1 * m . 

After simplifying, eq (5) becomes 

2 

r g±r—p sin 6 cos (<£— </>')+§- (6) 

This approximation of r gives a simpler form than that given by Silver [8] and also by Hansen 
and Bailin [3]. It should be noted that this particular form makes the following relatively 
simple analysis possible. 

Using the approximations (3) and (6), formula (2) gives: 

i~r/ N w e~ JkT C 2V C° n/ /N /fcfpsin* COS (*-*')-£l , - . 

EtrM^jt^J ) o *W')« L Hpf/p^. (7) 

If it is further assumed that the aperture illumination F(p<t>') is of circular symmetry, i.e., 
F(p,<t>') can be written as F(p); then E{rfi,<t>) will be independent of and can be denoted by 
E(r,6). Hence we have 



E(r,6)=jk — - P F(p)J (k sin P )« * W p ,/ p . 

P Jo 



(8) 



If a normalized variable 



H < 9 > 

is defined, then eq (8) can be written as 

E(r,B) =jka* (^Pj £ F(f) J («f) « ~' * " f#, (1 0) 

where 

u=ka sin (11) 

and 

to 2 



T= T 



(12) 



From eqs Ul) and (12) it is clear that u=y implies r sin d=a; therefore u=y means that the 
field point is at a distance, equal to the radius of the aperture, from the aperture-axis. 
It is more convenient for the following discussions to rewrite eq (10) as: 

E(r,B) = e-* r (ye j *-ir} f F(f) JoNV^ 1 "^ f#. ( 13 ) 
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3. Uniform Aperture Illumination 

For the case of uniform illumination, F(£) = l, eq (13) reduces to 

E(r, d) = e~^ (y* J ^)f Q Mufie' ?*-**> • fdf. (14) 

In the above equation, the integral may be denoted by Wq(j,u), 

W° (y,<u) =£ Mufie'i (1 " f2) m- (15) 

This integral has an imaginary part as well as a real part. It is closely related to the Lommel's 
functions of two variables [9]. The LommeFs functions of two variables can be expressed in 
integral forms as : 

U 1 (y,u)=y f 1 J M") cos J (1-f 2 ) ■ fif, 
Jo * 

U 2 (y,u)=y£ J («f) sin \ (1-f*) • tff. (16) 

Ui(y,u) and U 2 (y,u) are of order one and order two, respectively. 

By integrating by parts and using the recurrence relations of Bessel functions, it can be 
proved that Ui(y,u) and U 2 (y,u) m&y be expanded into the following two series: 

Ui ^ u)= (£) j ^ u) ~(£J Mu)+ (i) 5 Ji{u) ■■■' 

U 2 (y,u) =Q 2 J 2 (u) -g) 4 J*(u)+(l)° Mu) .... (17) 

Generalizations of these functions to higher orders are possible, but these generalizations will 
not be used in the present paper. 

In terms of Ui(y,u) and U 2 (y,u), we have 

my>u) jiihri +j ms.. (is) 

Using the above relation for Wo(y,u), formula (14) can be written simply as: 

E(r,d) = e- jkr \ye jW ~^) -W° (y,u). (19) 

This result is essentially the same as the first term in the series solution of Hansen and Bailin [3]. 
It can be seen from their figure 6 that the two results are close to each other within the region 
considered in this paper. 

4. Nonuniform Aperture Illumination of the Form (1— f 2 ) w 

If the aperture illumination is assumed to be of the form (1 — f 2 )", the field expression can 
be obtained by substituting F(£) = (l— f 2 ) w into eq (13), that gives 

E(r,6) = e-^ (ye'^f (W 2 ) *Jo(ut)e'* (W2) fdf. (20) 
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Now if we define 
then 



w n o(y,v>)=£ a-n n Mufie jl * (1 - 2) sdt (21) 

E(rfi)=e-** (ie^).Wl(y,u). (22) 



WJien n=^0, eqs (20) and (22) reduce to eqs (14) and (19), respectively. In other words 
eqs (20) and (22) include uniform illumination as a special case. For 717^0, these give the 
field distributions for the special tapered illumination (1 — f 2 ) n . The integration of the integral 
(21), or the evaluation of the function Wq (y,u), can be carried out by noting the following 
simple property: If differentiation of W$ (y,u) with respect to y is carried out, we have 

£wi(y,u)=%f\l-1*)'» 1 Mvt)fi <l - n tdt. (23) 

Equation (23) gives the recurrence relation, 

WF l {y,u)=~Wl(%v). (24) 

In terms of W° (y,u), we have 

W» (y,u)=Qy ^Wl(y,u) (25) 

or 

^<r,»>-0)"[£^«&^} (26, 

It is possible to derive expressions of W^{y,u) in terms of Lommel's functions of orders 
higher than two. But as far as numerical evaluation is concerned, it was found that the 
simplest way is to use eq (26) and the series expansions (17) for Ui(y,u) and U 2 (y,u). The 
series expansions of (<b n lby n )(Ui(y,u)ly) and (d n /dy n )(U 2 (y,u)/y), for n=0,l,2,3, and 4, are 
obtained by differentiating the power series of Uj(y,u)/y and U 2 (y,u)/y, term by term, n times 
with respect to y. 

5. Results 

In the last section, for the aperture illumination of the form (1 — f 2 ) n , the field expression 
is given by 



.7T— 7 



J ■ W" (y,u), 



Etr,e)=e-*\ye 2 ) ■ Wtfy.u), (27) 

where 

ha 2 
y=&- (28) 

and 

u=ka sin 0. (29) 

E(r,6) in eq (27) can always be rewritten in the following form: 

E(rfi)=Ae j2ir + (30) 

where A is the amplitude and 2ircj) is the phase angle in radians. 

Plots of the aperture illumination (1 — f 2 ) n for n = 0,l,2,3,4, are given in figure 3. The 
field distribution curves, both amplitude A and phase 27rc/>, are computed for the following 
cases (see fig. 4) : 
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Figure 3. Aperture illumination distributions, uniform and 
tapered. 




Aperture 



Figure 4. Fresnel region, in which Held distributions have 
been computed. 



(i) Field distributions along the circular arcs : 

1/T> 2 \ 1/D 2 \ 1/L> 2 \ /L> 2 \ ;/£> 2 \ 

where D—2a. The corresponding values for 7 are: 

3tt 7T 7T 

y=27T,— ; 7T, ->-, 0. 

It should be noted that the case r = ™ or 7=0 is the conventional far field result. 
For the convenience of comparison, the field distribution curves are normalized in the 
following way : Let 

E(r,0)=A o e* 2r +° (31) 



then the normalized E(r,d) for a fixed r is given by 



E(rfl) A e ~~^ n6 



(32) 



The normalized amplitude and phase distributions along the arcs are given in figure 5 to 

figure 14. 

(ii) Field distributions along the aperture-axis: 

In this and the following case, the distributions are not normalized but are plotted with 
the factor e~ jkr excluded. The two distribution curves, amplitude and phase, in the present 
case are useful in correlating the normalized distributions given in case (i). 
(iii) Field distributions along a line parallel to the aperture-axis and at a distance equal to 
the radius of the aperture : 
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Figure 5. Normalized amplitude distribution curves with 
aperture illumination— (1 — £* 2 ) , (constant). 
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Figure 6. Normalized phase distribution curves with aperture 
illumination = ( 1 — £ 2 ) °, (constant) . 




Figure 7. Normalized amplitude distribution curves with 
aperture illu m ination = (1 — £" 2 ) . 
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Figure 8. Normalized phase distribution curves with aperture 
illumination = (1 — £ 2 ) . 
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Figure 9. Normalized amplitude distribution curves with 
aperture illumination = (1 — f 2 ) 2 . 
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Figure 10. Normalized phase distribution curves with aperture 
illumination = ( 1 — £* 2 ) 2 . 
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Figure 11. Normalized amplitude distribution curves with 
aperture illumination = (1 — £ 2 ) 3 . 
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Figure 12. Normalized phase distribution curves with aperture 
illumination = (1 — f 2 ) 3 . 



Distribution curves for case (ii) and cases (iii) are given in figure 15 to figure 18. From 
these field distribution curves, the following properties are observed. 

(1) For the normalized field distribution curves along circular arcs: 

(a) For each aperture illumination of the form (1 — £ 2 ) n , the main beam width is essentially 
a constant in the range of r from }i(2D 2 /\) up to infinity. (Note: For a fixed direction, the 
parameter u is a constant.) 

(b) There is no perceptible difference in the normalized amplitude distribution curves at 
2D 2 /\ (conventional far field limit) and at infinity for the same aperture illumination, except 
near those directions where the infinity case has zero values or nulls. On the other hand, 
there are differences between the phase distribution curves; these curves change from a stepped 
form at infinity to a smooth continuous form at 2D 2 /\. 
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Figure 13. *' Normalized amplitude distribution curves with 
aperture illumination = (1 — £ 2 ) 4 . 
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Figure 14. Normalized phase distribution curves with aperture 
illumination = (1 — ^" 2 ) 4 . 



(c) The "side lobe level" increases as the distance r is reduced, but the amplitude distribution 
curve as a whole becomes "smoother." It is also true that the phase distribution curve 
becomes smoother too as r is reduced. 

(d) When the distribution curves of different aperture illuminations are compared, it is seen 
that the main lobe beamwidth increases as n is increased. At the same time, within a fixed 
measure of the angle 0, the number of side lobes is reduced. 

(e) The side lobe level is definitely reduced as n is increased. The ratio of the normalized 
maximum amplitudes of the first side lobes with n=A and n=0 is approximately 1 to 10. 

(f) As n increases, the difference in the normalized amplitude distribution curves at different 
r becomes smaller. 

(2) For the field distribution curves along the aperture-axis: 

(a) Along the aperture-axis, the amplitude distributions vary closely as 1/r in the range 
r=y 2 (2D 2 /\) to r= infinity for all cases studied. As for the particular case n=4, this is true 
even when r is reduced to %(2D 2 /\). 

(b) Along the aperture-axis, all the phase distributions, excluding the exponential factor 
e~ jkr , are approximately straight lines, i.e., the phase varies linearly with respect to 1/r. 

(3) For the field distribution curves along a line parallel to the aperture-axis and at a 
distance equal to the radius of the aperture; 

(a) In the range, r=2D 2 /\ to r=infinity, both the amplitude and phase distributions have 
approximately the same values as those obtained along the aperture-axis. 

(b) At r=%(2D 2 /\), the ratio of the corresponding amplitudes for the present case and case 
(ii) is roughly one to four for all cases computed. 

6. Generalization to Nonuniform Aperture Illuminations of Other Forms 

In the last section, only a special type of nonuniform illumination of the form (1 — f) n 
was given. It is possible to generalize the analysis to include nonuniform illuminations of 
olhcr forms either through expansions by the Schmidt orthogonalization process [5] or by- the 
following simple algebraic procedure. For any circularly symmetric aperture illumination, it 
is always possible to approximate it by an even polynomial of the variable f. Such poly- 
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Figure 15. Amplitude distribution curves along the aperture 
axis with aperture illumination=(l — £ 2 ) n . 
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Figure 16. P/iase distribution curves {excluding the factor e" jkr ) 
aZongf 2/ie aperture axis with aperture illumination = (1 — £ 2 ) w . , 
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Figure 17. Amplitude distribution curves along a line parallel 
to the aperture axis and at a distance a from the aperture 
axis with aperture illumination = (1 — f 2 ) n . 
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Figure 18. Phase distribution curves (excluding the factor e* jkr ) 
along a line parallel to and at a distance a from the aperture 
axis with aperture illumination=(l — £ 2 ) n . 
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nomials are of the following form: 

a„{ 2n +a n ^ 2(n - l) + • • • +aif 2 +«o. (33) 

It is clear that such a polynomial can always be transformed into the following form: 

b n (l-n a +K-i(l-n n - 1 + • • • +bi(l-n+K (34) 

For a given set of a's, the corresponding 6's can be uniquely determined by solving a system of 
linear simultaneous algebraic equations. In particular, for the case of ft=4, the coefficients 
&'s obtained in terms of the a's are as follows: 

b±=a 4 

bz=— a- 3 — 4a 4 

62=^2+3^3+6^4 

b l =—a i — 2a> — 3a 3 — 4a 4 

b =a + ai+a 2 +a 3 +a 4 . (35) 

When 6 , 6 ]? . . . , 6 4 are determined as above, the field expression is then given by 

E(r, e)=e~ jkr \yJ*) [ * 2 6,0 -f 2 ) n Mufie* {1 ~^^ (36) 

Jo //=() 

In terms of IVd' (\,w), (36) gives 

E(r, e)=e-*» (ye 1 **) ± b n W?(y,u) (37) 

Equation (37) can then be evaluated by using the results already obtained in this paper. 
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